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ABSTRACT 

A  five-level  baroclinic  model  for  numerical  prediction  of  height 
fields  at  1000,  850,  700,  500,  and  300  mb  is  developed  and  programmed. 
Several  forms  of  the  model  are  tested  to  determine  the  contribution 
of  the  various  terms  of  the  vorticity  equation.   Forecasts  from  three 
different  methods  are  compared.   The  model  shows  utility  in  prognosis. 
Although  there  are  individual  differences  among  the  methods,  no  one  of 
the  three  is  significantly  better  over  the  entire  hemisphere  on  the 
data  tested  here. 
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Table  of  Symbols 

Q  the  upward  component  of  the  apparent  gravitational  acceleration 

D  the  altimeter  correction    P~  *?  ~Zy 

Z  the  height  of  an  isobaric  surface 

Zp  the  pressure  altitude  in  the  standard  atmosphere 

-f  the  Coriolis  parameter,  Ofl   siv-s^   where  (p    is  the  geograph- 
ical latitude 

T  the  mean  value  of  the  Coriolis  parameter,  2.   il  Sir>  H  5 

$  the  geopotential 

£  the  relative  vorticity   £  =  -|-v  -  ^ 

A  the  absolute  vorticity   >?;  X     ■+  f 

~\\f  the  stream  function  for  the  non-divergent  component  of  velocity 

X  the  velocity  potential  for  the  divergent  component  of  velocity 

CO  the  vertical  velocity  of  air  in  pressure  coordinates 

P  the  atmospheric  pressure  in  millibars 

0"  the  static  stability     G"  ■  -  q"  ^p 

K  the  universal  gas  constant 

d  grid  distance,  381  km  at  60N 

W  map  scale  factor 

V  del  operator  on  a  constant-pressure  surface 

V  horizontal  Laplacian  operator  on  a  constant-pressure  surface 
J  horizontal  Jacobian  operator   J  (A,  B)=  "^  yr;  ~^y  y7" 

V  finite  difference  Laplacian 

JJ  finite  difference  Jacobian 

J(A3)=[(A(+,J-A,uXB,J-BLj-y(Ap,1-/\)>,yBM,j-I3(,(J 
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1.   Introduction. 

Until  very  recently,  the  barotropic  models  have  been  the  principal 
source  of  numerical  weather  prognoses.   Their  shortcomings  in  areas  of 
baroclinic  development,  plus  the  advent  of  sophisticated  digital  computer 
systems,  have  spurred  interest  in  baroclinic  models.   It  is  the  purpose 
of  this  research  to  develop  and  program  a  five-level  baroclinic  model  and 
to  investigate  the  contribution  of  the  various  terms  of  the  vorticity 
equation  to  the  prognosis  of  the  pressure  field. 


2.   Background. 

In  the  model  developed  here,  prognosis  is  accomplished  by  means  of 
the  vorticity  equation,  which,  with  pressure  as  the  vertical  coordinate, 
can  be  written 

Basic  to  the  derivation  of  this  equation  is  the  Helmholtz  theorem  which 
expresses  the  wind  as  the  sum  of  a  divergent  and  a  non-divergent  compon- 
ent: 

.  V-Ik  x  Vf  +  V*-.  (2) 

The  potential  function  X.  for  the  divergent  part  of  the  wind  is  obtained 
through  solution  of  the  continuity  equation 


V  X  + 


~duo 
<>p 


=  o 


(3) 


The  stream  function  ~y    for  the  non-divergent  portion  is  approximated  by 
solution  of  the  balance  equation 

as  described  by  G.  'Arnason  [lj,  where  d)  is  the  geopotential.   Use  of 
the  stream  function  in  the  vorticity  equation  enables  one  to  avoid  the 
geos trophic  approximation  and  the  concomitant  spurious  anticyclogenesis. 
The  vertical  component  of  relative  vorticity,  C  ,  is  the  Laplacian 
of  the  stream  function 


C-7T 


(5) 


and  absolute  vorticity,  '[   ,  is  the  sum  of  this  and  the  vertical  component 
of  the  earth's  vorticity,  f. 

v  =  r  +  f  (6) 

Vertical  velocity,  60,  is  computed  by  the  diagnostic  equation 


v-^^T  !^--#kJUT)-j(z,^T)4j(-u) 


(tr  apl  f 


(7) 


developed  by  Clarke  and  Lawniczak  [3]  .   The  forcing  function  is  calculated 
with  the  geostrophic  approximation.   This  results  in  some  inconsistencies 
in  the  model  but  to  avoid  reprogramming  the  entire  equation,  we  have  gone 
ahead  with  comparison- type  experiments.   Vertical  velocity  is  assumed  zero 
at  100  mb,  but  other  boundary  conditions  are  among  the  options  available 
to  the  user  of  this  computer  program  as  follows: 

a)  the  lower  boundary  can  be  level  or  include  terrain; 

b)  vertical  velocity  can  be  set  equal  to  zero  at  the  lower  boundary, 
or  computed  with  or  without  terrain  effects  and  with  or  without  the 
effect  of  friction; 

c)  the  stability  parameter  (T  can  be  a  point  variable,  a  function  of 
pressure  alone,  or  a  constant;  and 

d)  the  maximum  acceptable  residual  in  the  3-dimensional  relaxation 
can  be  selected  by  the  user. 

In  the  vorticity  equation,  the  contributions  from  the  twisting  term, 

/  loo  ??t"       a«o  yyr  \ 

( "a-*-  w  3p  ~~  Yy    by  3p  J   an<*  fc^e  vertical  advection  of  vorticity  term, 
U)  -£-*-    are  usually  considered  small,  and  some  investigations  have 
shown  that  they  tend  to  compensate  each  other  [2],   Also,  integral  theo- 
rems show  [4j  that  to  avoid  the  mean  generation  of  vorticity  in  a  closed 


volume,  the  two  terms  should  be  included  or  neglected  together.   Both 
have  been  omitted  from  the  initial  investigation  reported  here. 

A  further  simplification  involves  omission  of  the  advection  of  abso- 
lute vorticity  by  the  divergent  portion  of  the  wind,  5—  sr~ *-  +  -r- -   -.-^  . 
For  consistency,  it  is  necessary  then  also  to  substitute  the  mean  value 
of  the  Coriolis  parameter,  -f  ,  for  the  absolute  vorticity  in  the  diver- 
gence term,  f  -|^-  . 

In  the  equations  that  follow,  the  stream  function  has  been  redefined 
so  as  to  be  in  units  of  height, 


r-fE. 


(8) 


The  resulting  forms  of  the  vorticity  equation,  expressed  in  finite- 
difference  notation,  are 

TT^m     7^pJJ(t.^J  +  f  apt  <9> 


We>  w^l  __  -fl/Vyf.Tr/r    7>\     _^/^^4^^Uft^        CIO* 
and   the  complete  equation 

where    cj   is   the  grid  distance  and  /W\  is   the  ratio  of  this   distance   to   the 
corresponding  geographical   distance. 


Horizontal  spatial  operators  are  approximated  by  the  usual  centered 
differences.   To  approximate  the  vertical  derivatives  in  the  vorticity 
equation  a  Lagrangian  parabola  was  fitted  through  three  points.   For  a 
centered  difference  the  equation  is 


,   CP.-PQ    A  ,  UP.-P.-PJ  A  ,    (P.-P.)    A       Q2, 


^P  p-p,  CP.-M(P.-ft)    (PrOCPrPO    (P,  P.)(P,-P,)' 

All  computations  were  done  on  the  Control  Data  Corporation  model 
1604  computer,  using  the  U.  S.  Navy  Fleet  Numerical  Weather  Facility 
octagonal  grid.   This  grid,  centered  at  the  pole  and  inscribed  in  the  ION 
latitude  circle,  consists  of  1977  grid  points  with  a  grid  size  of  381  km 
at  60N. 

The  week  of  22  November  1962  was  a  period  of  considerable  cyclogene- 
sis  in  dense-data  areas,  and  therefore  excellent  for  testing  a  baroclinic 
model.   Initial  data  used  were  the  U.  S.  Navy  Fleet  Numerical  Weather 
Facility  height  analyses  at  1000,  850,  700,  500,  and  300  mb  for  1200Z  on 
22  November  1962.   The  height  analyses  for  0000Z  and  1200Z  23  November 
were  used  in  the  evaluation  of  the  predictions. 


3.   Procedures. 

In  the  diagnostic  60 -equation  for  the  initial  tests  discussed  here, 
the  boundary  conditions  used  are: 

a)  vertical  velocity  equal  to  zero  at  100  mb; 

b)  vertical  velocity  equal  to  zero  at  the  lower  boundary; 

c)  no  terrain  effects  at  the  lower  boundary;  and 

d)  the  stability  parameter,  g-  ,  is  a  function  of  pressure.   The  pro- 
gram is  used  in  "forecast  mode";  that  is,  temperatures  are  computed  hydro- 
statically,  rather  than  taken  from  objective  analysis.   The  maximum 
residual  in  the  3-dimensional  relaxation  is  1  x  10   mb  per  second. 

The  sequence  of  calculations  programmed  can  best  be  explained  with 
reference  to  figure  1.   Input  data  in  every  case  are  the  U.  S.  Navy  Fleet 
Numerical  Weather  Facility  height  analyses  at  the  five  standard  levels, 
modified  to  ensure  vertical  consistency. 

In  the  first  method  with  equation  (9),  stream  functions  are  used 
throughout.   The  inversion  to  height  fields  [l]  ,  effectively  the  reverse 
of  the  balance  equation,  is  thus  bypassed  and  stream  functions  are  used 
for  input  data  for  both  diagnostic  and  prognostic  equation.   This  is 
hereafter  referred  to  as  method  I.   The  resulting  24-hour  prognoses  are 
displayed  on  charts  6  through  10  and  discussed  in  section  4  below. 

Since  the  geostrophic  approximation  is  used  in  the  forcing  function 
of  the  CO -equation,  and  temperatures  are  computed  hydrostatically ,  there 
might  be  some  reason  to  suspect  the  vertical  velocities  computed  from 
height  fields  could  be  more  accurate  than  those  from  stream  functions.   A 
completely  geostrophic  model  was  programmed,  using  pressure-height  data 
throughout.   This  method  showed  the  expected  spurious  anticyclogenesis 
even  at  12  hours  and  was  carried  no  further. 
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The  complete  process  diagrammed  in  figure  1  has  been  used  in  the  re- 
mainder of  this  investigation.   In  the  first  time  step  the  initial  data 
(five  vertically-consistent  height  fields)  are  used  as  input  to  the  bal- 
ance equation  and  to  the  00  -equation.   From  the  former  are  obtained  stream 
functions  for  the  five  levels;   the  latter  yields  vertical  velocities  at 
the  same  levels.   With  these  stream  functions  and  vertical  velocities  the 
vorticity  equation  (equation  (9)  for  method  II,  equation  (10)  for  method 

III)  gives  — -  ,  the  local  change  in  stream  function,  at  the  five  levels . 

At 

When  these  are  added  to  the  initial  values,  the  first  forecast  time  step 
is  complete.   As  shown  in  figure  1,  these  new  E  values  are  now  input 
fields  for  both  the  vorticity  equation  and  the  balance  equation;   the 
latter  being  used  for  inversion  to  pressure  height  as  required  by  the 
6J -equation. 

Subsequent  time  steps  follow  the  same  routine,  except  in  the  applica- 
tion of  ^E  ,  which  is  done  according  to  the  centered-dif ference  equation 
At 

£„.,  +i(-Hl,^  =EK41  .  (13) 

In  order  to  minimize  the  errors  inherent  in  the  initial  forward 
time  step,  the  first  time  increment  is  only  15  minutes.   Centered  time 
steps  are  then  used  to  complete  successive  15-minute  prognoses  to  one 
hour,  and  finally  1-hour  prognoses  to  12  hours.   The  fields  are  then 
smoothed  lightly  to  remove  any  small  perturbations  resulting  from  trunca- 
tion or  round-off  errors,  and  the  process  is  repeated  to  24  hours. 


4.    Results. 

In  review:  method  I  uses  stream  functions  throughout,  and  with 
equation  (9)  as  the  prognostic  equation  excludes  the  advection  of  abso- 
lute vorticity  by  the  divergent  portion  of  the  wind  and  substitutes  f 
for  7l    in  the  divergence  term;  method  II  inverts  the  stream- function 
fields  to  pressure  heights  at  every  time  step,  but  retains  the  other 
simplifications  of  method  I;  method  III  also  inverts  to  pressure  height 
at  each  time  step,  but  uses  the  more  complete  equation  (10)  as  the  prog- 
nostic equation. 

Charts  1-5  are  the  pressure  height  analyses  at  the  beginning  of  the 
period.   The  next  15  charts  are  the  results  of  the  three  24-hour  fore- 
casts: method  I  in  charts  6-10;  method  II,  charts  11-15;  and  method  III, 
charts  16-20.   Charts  21-25  are  the  observed  height  analyses  at  the  end 
of  the  forecast  period.   Chart  26  is  the  U.  S.  Navy  Fleet  Numerical 
Weather  Facility  barotropic  forecast  at  500  mb  for  the  same  period „ 
Table  1  shows  the  contour  origin  and  interval  for  all  the  charts. 

In  general,  corresponding  charts  from  the  three  methods  are  quite 
similar.   They  all  retain  a  typical  meteorological  appearance  and  the 
various  systems  show  a  reasonable  vertical  structure.   Further  indica- 
tion of  their  similarity  is  contained  in  table  2.   The  "pillow"  referred 
to  is 

Pillow   =  ""  (A  "  &)"  (16) 

M77 

and  the  RMSE  is  the  root-mean-square  error  calculated  in  the  usual 
fashion: 
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y    [(A-Ej-p,u-w]m 

RMSE   =  |  *ZH .       (17) 


1^77 

Nevertheless,  a  level-by- level  comparison  of  the  three  forecast 
methods  shows  some  differences  among  them.   At  1000  mb,  for  instance,  the 
two  low  centers  over  North  America  merged  and  deepened,  while  the  trough 
moved  nearly  20  east.   All  methods  forecast  the  merging  and  consider- 
able trough  movement,  with  method  I  giving  the  best  placement  of  the 
trough.   During  the  forecast  period  a  low  center  developed  in  the  area 
of  weak  gradient  southwest  of  the  Kamchatka  peninsula.   This  development 
was  forecast  by  all  methods,  although  with  too  much  cycloganesis. 
Method  I  was  again  the  best.   In  no  case  is  there  any  evidence  of  spuri- 
ous anticyclogenesis,  even  at  low  latitudes. 

The  forecasts  at  850  and  700  mb  show  movement  consistent  with  that 
forecast  at  1000  mb.   By  700  mb  the  tilt  of  the  baroclinic  systems  has 
become  apparent.   A  good  example  is  the  stationary  low  center  over  the 
Northwest  Territories,  which  deepened  by  26  meters.   All  forecasts  held 
it  stationary,  but  methods  I  and  III  forecast  two  and  seven  meters  deep- 
ening, respectively,  while  method  II  forecast  a  2-meter  filling.   The 
weakening  of  the  low  center  east  of  Greenland  was  forecast  by  all  methods, 
but  best  by  method  II. 

At  500  mb,  since  the  chart  of  the  U.  S.  Navy  Fleet  Numerical  Weather 
Facility  barotropic  forecast  has  been  included,  there  are  four  charts  to 
compare.   There  is  a  striking  similarity  between  the  barotropic  progno- 
sis and  that  from  method  I.   The  method  II  and  III  prognoses,  themselves 
much  alike,  are  quite  different  from  the  other  two.  Method  I  and  the 
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USNFNWF  model  both  forecast  the  movement  and  reorientation  of  the  trough 
over  Eastern  United  States,  and  the  weakening  and  displacement  of  the 
low  center  near  Hudson  Bay.  Methods  II  and  III  show  the  trough  lagging 
far  behind  its  true  position,  and  retain  the  low  as  a  stationary  system. 
Method  III  even  deepened  it  slightly.   In  contrast,  methods  II  and  III 
produced  much  better  forecasts  than  the  others  in  the  western  Europe- 
Mediterranean  area. 

At  300  mb  all  three  methods  agree  in  general  on  the  movement  and 
development  of  the  major  systems.   This  development  is  in  some  cases  out 
of  phase  with  the  actual  development,  but  no  one  method  can  be  selected 
as  better  than  the  rest.   The  initial  and  verifying  charts  at  300  mb 
show  more  effect  from  the  process  of  making  the  data  vertically  consist- 
ent than  was  evident  at  lower  levels. 

Referring  again  to  table  2,  there  seems  to  be  no  best  method  among 
the  three  tested.   Certainly  the  addition  of  the  advection  of  absolute 
vorticity  by  the  divergent  wind  and  the  use  of  absolute  vorticity  in  the 
divergence  term  produced  little  change  and  no  over-all  improvement  in 
this  particular  forecast.   The  effect  of  the  inclusion  of  the  vertical 
advection  of  absolute  vorticity  and  the  twisting  term  remains  to  be 
tested. 
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Table  1 


Level 
(mb) 

1000 

850 

700 

500 

300 


Contour  or 
(meters 

i$in 
) 

Cont 

:our  interval 
(meters) 

120 

60 

1440 

60 

3060 

60 

5580 

60 

9120 

120 
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Table  2 


Me  thod 


Level 


Pillow 


RMSE 


(mb) 

(meters) 

(meters) 

I 

1000 

-  7.9 

43.6 

850 

-  1.8 

37.6 

700 

-  1.8 

41.8 

500 

-  1.2 

51.3 

300 

13.4 

80.3 

II 

1000 

-11.3 

43.3 

850 

-  2.7 

40.8 

700 

0 

42.7 

500 

-   .6 

53.4 

300 

11.9 

77.3 

III 

1000 

-  9.1 

49.1 

850 

-  3.4 

43.0 

700 

-  2.4 

43.6 

500 

-  1.5 

54.4 

300 

13.4 

80.0 

USNFNWF 

500 

-  4.9 

47.0 

Persistence 

500 

-  3.3 

57.4 
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APPENDIX  I 


The  following  charts  were  computed  during  the  course  of  this  study, 
They  are  included  here  for  easy  reference. 


Charts  1-5 
Charts  6-10 
Charts  11-15 
Charts  16-20 
Charts  21-25 
Chart  26 


Analyses,  12Z  22  Nov  1962 
24-hour  prognoses,  method  I 
24-hour  prognoses,  method  II 
24-hour  prognoses,  method  III 
Analyses,  12Z  23  Nov  1962 
24-hour  prognosis,  USNFNWF 
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■       •  _  SCALE:  1: 60,000,000 V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


Chart  N0..6-B 


858    HT       PMPL     "12Z..22  NOU-62    CHART    2 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60°  TMORTH  LATITUDE 
' SCALE:  1:60,000,000        V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


CHART  NO    6-B 


700    HT       AMflL    /'12Z--22  NOU  G2  CHART  3 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60°  'NORTH  LATITUDE 

•_ ■     '  ■        ; scale: _U«o.oqo!ooq__  v. 


FLEET  NUMERICAL  WEATHER  FACILITY 
-JL       MONTEREY,  CALIFORNIA 


300     HT       PINAL      .  12Z-  22  NOU  €2      CHARTS 

67*7       * 
PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60"  ><ORTH  LATITUDE 
j. ___ SCALE:  1:60,000,000 y. 


1000  HT  21  HR  PROG  FROM  -  12Z  22 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60"  NORTH  LATITUDE 
SCALE:  1;  60,000.000  V 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


850     HT  21  HR  PROG  FROM     12Z  22  NCJiU  62   CHART  7 


■3z*&      * 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60°  NORTH  LATITUDE 
SCALE:   1:60,000.000  if. 


FLEET  NUMERICAL  WEATHER   FACILITY 
MONTEREY,  CALIFORNIA 


700  HT  24  HR  PROG  FROM  12Z  22  NQU  G2  CHART  8 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  50°  IMORTH   LATITUDE 
SCALE.   1:60.000,000  V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY.  CALIFORNIA 


CH«RT  NO    6  B 


500     HT  24  HR  PROG  FBOM  ■  12Z  22  N^U- G2    CHART    9 


PROJECTION:  POUR  STEREOGRAPHIC-TRUE  AT  60"  »»ORTH  LATITUDE 
SCALE:  1:60.000000  _.__  V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


308     HT  24  HR  PROG  FROM     12Z  22  NQU  62   CHART    10 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60°  NORTH  LATITUDE 
SCALE:    1:60.000.000  V 


f  FLEET  NUMERICAL  WEATHER  FACILITY 

S  MONTEREY.  CALIFORNIA 


CHART  NO.  6-B 


*  .       I    •  , 


?m     HT  21  HR  PROG  FROM     12Z  22NQU.G2     CHART    13 


;nH  "       ■:■:■■: 

PROJECTION:  POLAR  STEHEOORAPHIC-TRUE  AT  60' "NORTH  LATITUDE 
'  ^CALE:  1:60.000.000  V- 


S&d    HT  24:  HR  PROG  fBOM  "  12Z  22  Ni 

•  '  ' .  737*>3  + 

PROJECTION:  POLAR  STEREOORAPMIC-TRUE  AT  «T  >IORTH  LATITUDE 
__.: ■       '    ■  SCALE:  1:60,000,000  V- 


300     HT  24  HR  PROG  FBOM  '  12Z  22  NdiU- G2    CHART   15    * 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60'  *IORTH  LATITUDE 
9CALE:  1:  $0,000,000  if. 


iee&'-.KT  2«±.HR  PROG  FROM  "  12Z  22  .NQU- 62    CHART   16 


POLAR  STEREOGRAPMIC-TRUE  AT  W  NORTH  LATITUDE 
SCALE:  1:60.000.000  V 


ew  NP..M 


7®d     HT  21  HR  PROG  FROM     12Z  22~N0U  G2    CHART   18 


3i_  T*S 
PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60'  NORTH  LATITUDE 
SCALE:    1:60,000,000  V. 


f  FLEET  NUMERICAL  WEATHER   FACILITY 

MONTEREY,  CALIFORNIA 


CHART  NO.  6-B 


300     HT  21  HR  PROG  FROM  ■  12Z  22  MOU  62   CHART    20 


32  #}      * 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60'  NORTH  LATITUDE 
SCALE:   1:60.000.000  V 


Chart  No  6-B 


100&    HT       ANAL     *  122:23  WU  62     CHART   21 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60"  NORTH  LATITUDE 
SCALE:   1:60,000,000  V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


50&     HT       ANAL     "  12Z * 23  NOU  <G2     CHART  24 


&>*'    *        ■.'■'.* 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60'  KORTH  LATITUDE 
SCALE:   1:60,000,000 if. 


.^•* 


JtiHfl  985 


iJAOS 


30&     HT       PMflL     *  12Zj 23  N0U-G2    CHART    2  5 


C7*Z      * 

PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60" 'NORTH  LATITUDE 
SCALE:   1:60,000,000 V- 


FLEET  NUMERICAL  WEATHER  FACILITY 
MONTEREY,  CALIFORNIA 


*  •  .    '    ' » 


'50&     HT  24  HR  PROG  FROM  "  12Z  22  -N 


PROJECTION:  POLAR  STEREOGRAPHIC-TRUE  AT  60"  NORTH  LATITUDE 
SCALE:   1:60,000,000  V- 


